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(N Abstract 

The concept of the broadband epsilon-near-zero meta-atom consisting of layered stacks with 

specified metallic filling ratio and thickness is proposed based on the Bergman spectral representa- 

tion of the effective permittivity. The step-like metal-dielectric multilayer structures are designed 

c/3 to achieve realistic broadband epsilon-near-zero meta-atoms in optical frequency range. These 
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meta-atoms can be integrated as building blocks for unconventional optical components with ex- 
otic electromagnetic properties over a wide frequency range, such as the demonstrated broadband 
directional emission and phase front shaping. 
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I. INTRODUCTION 

Metamaterials are artificially structured materials with subwavelength components, which 
can be designed to create extraordinary macroscopic electromagnetic properties that do not 
exist in nature 1-11 . Among all kind of metamaterials, the metamaterials with near- zero 
permittivity (epsilon-near-zero, ENZ) emerge into the focus of the extensive exploration 
due to their anomalous electromagnetic features at microwave and optical frequencies 12-23 . 
However, all ENZ metamaterials suffer from single frequency response for the near-zero 
permittivity, which gives a huge disadvantage for practical applications. 

In this paper, we further develop the theoretical strategy 24 for mathematically designing 
the broadband ENZ metamaterials based on the Bergman spectral representation 25,26 of the 
effective permittivity, and propose the realistic broadband ENZ meta-atoms by applying 
the step-like metal-dielectric multilayer structures. In quasi-static conditions, the Bergman 
spectral representation expresses the effective permittivity of a composite medium in terms 
of a series of residue-singularity couples. Therefore, the broadband ENZ metamaterial can 
be flexibly designed in a dimensionless spectral space, while the geometric information of 
the metamaterial, including the dimensions and the filling ratio of each component, can 
be obtained through a proper inverse problem. According to this strategy, the step-like 
metal-dielectric multilayer structures are proposed to realize the broadband ENZ meta- 
atoms in optical frequency range. Functional optical devices constructed by the realistic 
broadband ENZ meta-atoms are also demonstrated, and exotic electromagnetic properties 
of such devices are explored including the broadband directional emission and phase front 
shaping. 

II. DESIGN OF BROADBAND ENZ META-ATOM 

As depicted in Fig. 1(a), the proposed broadband ENZ meta-atom is a two dimensional 
three-layer stack with the dimensions of 100 nm x 25 nm. Each layer of the meta-atom is 
a metal-dielectric mixture with gold (Au) as the metallic inclusion and silica (SiC^) as the 
dielectric host medium. The permittivity of gold follows the simple Drude model Eau = £00 — 
Up/iuiuj + i^)) with dielectric constant 6^ = 5.7, plasma frequency u p = 1.3666 x 10 16 rad/s, 
and damping constant 7 = 3 x 4.0715 x 10 13 rad/s. Here the damping constant is of the 



value three times higher than the bulk value in order to account for surface scattering, grain 
boundary effects, and inhomogeneous broadening of the metal thin film. The permittivity 
of silica is £si0 2 = 2.1338. In quasi-static conditions, the Bergman spectral representation 
expresses the effective permittivity of the meta-atom as 



£e = e Si02 (l - F(s)) 
with respect to the variable s = £si0 2 /(£si0 2 — £Au)- The spectral function F(s) reads 
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which is analytic everywhere with positive residue Fi except for three simple singularities Sj 
confined as ^ S\ < s% < S3 < 1. Correspondingly, based on the geometric structure of the 
stack, including the filling ratio fi of the gold inclusion in each layer and the thickness di of 
each layer, the effective permittivity e x of the whole stack reads 
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with the total stack thickness d = Yli=i di- Through the effective permittivity, the connec- 
tion between the singularities and residues in the spectral space [Eq. (1)] and the geometric 
structure of the stack in the physical space [Eq. (3)] can be set up as 
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Clearly, the variable s is only a function of the frequency u as s = s(u), which implies 
the Bergman spectral representation on the left hand side of Eq. (4) describes the effective 
permittivity of the stack only as a function of the frequency u. Therefore, for designing the 
broadband ENZ meta-atom, it is simple to request that 
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with respect to a series of specified frequencies Uj (j = 1, 2, 3) in a wide frequency range 
[uj a , u>b\- Regarding the mathematical structure of Eqs. (3) and (5), the value of the singular- 
ity Si in Eq. (5) can be set as S\ = and Sj = Re(s(o;j)) for i — 2, 3 with respect to another 



series of specified frequencies uji within [cu a , u;&]. The specified frequency Ui and Uj can be 
simply arranged as an equally spaced array based on the following confinement 

U a = U j=1 < U i=2 < U j=2 < U i=3 < U j=3 = U b . (6) 

Therefore, the residue Fi can be fully determined by Eq. (5). With the value of the residue 
Fi and the singularity s,, the filling ratio fi and the thickness di of the broadband ENZ 
meta-atom can be achieved from Eq. (4). It is worth noting that this strategy is suitable 
for designing multilayered meta-atom for realizing other constant permittivity over a wide 
frequency range in principle, and the number of stack layers is not limited. In order to confine 
the optical loss within a moderate range, the operating frequency range of the broadband 
ENZ meta-atom in Fig. 1(a) is set between u a = 537.42 THz and w& = 589.62 THz, and the 
filling ratio fi of the gold inclusion and the thickness di of each layer are designed and listed 
in Table 1. According to the results, the realistic meta-atom, shown in the right panel of 

TABLE I. Filling ratio of the gold inclusion and the thickness of each layer 



Layer fi di (nm) 

1 0.216329 43.6884 

2 0.188523 19.2691 

3 0.163955 37.0425 

Fig. 1(a), is also designed as a three-layer stack with the same thickness di of each layer, but 
each layer is replaced by a gold-silica multilayer structure, where the height hi of the gold 
is hi = fi ■ h (i = 1,2,3) with respect to the same filling ratio fi of the gold inclusion. It 
is worth mentioning that in practical fabrication, the feature size of the fabricated thin film 
has the accuracy of nanometer, and a detailed analysis on the robustness of the designed 
parameters is presented in the Appendix. 

The effective permittivity e x of the effective medium meta-atom is plotted in Fig. 1(b) 
with the red curves standing for the real part and blue curves standing for the imaginary 
part (the optical loss). Moreover, the solid curves indicate the theoretical value of e x based 
on Eq. (5), while the dashed curves indicate the retrieved value, which is calculated from 
the scattering parameters (Sn and S21) of the effective medium meta-atom from the full- 
wave finite-difference time-domain (FDTD) method. Five equally spaced frequencies (ENZ 



frequency), at which the value of e x equals to zero, are listed at the right bottom of the 
figure and marked as yellow dots on the curve. ENZ frequencies i, iii, and v represent 
the frequency Uj, while ii and iv stands for the frequency Wj. It is clear that the FDTD 
retrieved results matches the theoretical results well, with only small fluctuations caused 
by the inevitable phase variation of the reflected electromagnetic wave from different stack 
layers in the effective medium meta-atom, which is not included in the theoretical strategy 
Additionally, the value of the effective permittivity e y of the effective medium meta-atom is 
e y ~ 2.8 + O.OOli in the operating frequency range. 

The upper panel of Fig. 2 gives the distributions of the absolute value of the local electric 
field inside the effective medium meta-atom with a height of 250 nm at the five ENZ frequen- 
cies with the electromagnetic wave directly illuminating from the bottom of the meta-atom, 
while the lower panel of Fig. 2 shows the electric field inside the realistic meta-atom array 
with the same height (including 10 realistic meta-atoms along the y-direction). As indicated 
in Fig. 2, the electric field is strongly confined and guided through the meta-atom along 
different paths with respect to different ENZ frequencies, due to the continuity of the elec- 
tric displacement with near- zero permittivity. The properly designed filling ratio of the gold 
inclusion in each stack layer ensures the zero effective permittivity around each specified 
ENZ frequency. Then the corresponding thickness of each stack layer modifies the inter- 
actions between the adjacent layers, resulting in the broadband ENZ response across the 
desired frequency range. Additionally, the similar distribution of the electric field inside the 
effective medium meta-atom and the realistic meta-atom array, although disturbed by the 
inevitable scattering caused by the realistic structures, informs that the realistic meta-atom 
also possesses a similar broadband ENZ response that is coincident with the theoretical 
design. 

III. ANALYSIS OF BROADBAND ENZ OPTICAL DEVICES 

Functional optical devices with broadband response corresponding to the broadband ENZ 
property can be constructed by the realistic broadband ENZ meta-atom. For example, 
Fig. 3(a) schematically gives a prism made of the meta-atoms in air, with an oblique up- 
per surface of angle 9 for realizing the directional emission of electromagnetic waves. The 
incoming electromagnetic wave with vertical wave vector fc in and Poynting vector S- in can 



be collimated inside the prism as it crosses the lower interface from air, while the refracted 
out-going electromagnetic wave with the wave vector fc out and Poynting vector S ont will be 
normal to the upper interface. Physically, the directional emission is resulted in the conser- 
vation of the tangential component of the wave vectors at the upper interface between the 
prism and air, which can be explained by the flat elliptical iso-frequency contour (IFC) of 
the meta-atom compared with the circular IFC of air with radius k$ = 2tt/X in Figs. 3(b) 
and 3(c). Clearly, the propagating direction of the out-going electromagnetic wave will be 
quite close to the normal direction (k' -direction) of the upper interface of the prism without 
the consideration of material loss [Fig. 3(b)], but slightly bend away when the material loss 
is taken into account [Fig. 3(c)]. 

Correspondingly, Fig. 4 gives the simulation results of the directional emission for prisms 
made of effective medium meta-atoms [Fig. 4(a)] and realistic meta-atoms [Fig. 4(b)] with an 
oblique interface of 9 = 15° at the designed ENZ frequency of 537.42 THz. The distribution of 
the electric field amplitude E x is displayed, with the power flow indicated by hollowed arrows. 
Strong reflection occurs at the bottom of the prism due to the mismatched impedance 
(analysis about the impedance matching to free space is given in the Appendix), while the 
near-zero phase variation of the propagating electromagnetic wave inside the prism and the 
small angle of refraction of the out-going electromagnetic wave at the upper interface of 
the prism are observed due to the near-zero permittivity. Furthermore, the variation of the 
angle of refraction with respect to different frequencies [Fig. 4(c)] reveals the broadband 
ENZ response of the realistic meta-atom made prism, through the comparison among the 
theoretical results obtained from the Snell's law (blue curve), the simulation results of the 
effective medium prism (blue dots), and the simulation results of the realistic meta-atom 
made prism (black circles). It is worth mentioning that the deviations of the simulation 
results of the realistic meta-atom made prism is due to the inevitable material loss and 
scatterings caused by the internal structures, and the average angle of refraction of the 
realistic meta-atom made prism is about 6.3° (red line) in the designed operating frequency 
range, representing an effective group index of n g = 0.42. 

Because of the relatively small phase variation caused by the near-zero permittivity, the 
realistic broadband ENZ meta-atom can also be used to build optical devices for shaping 
the phase front of electromagnetic waves as the S"-shape lens example shown in Fig. 5. 
The example is demonstrated at the ENZ frequency of 537.42 THz for S-shape lenses made 



of effective medium meta-atom [Fig. 5(a)] and realistic meta-atom [Fig. 5(b)], with the 
results indicated by the distribution of the electric field amplitude E x . Clearly, the uniform 
phase pattern of the electromagnetic wave inside the S"-shape lens due to the near-zero 
permittivity leads to a phase front of the out-going electromagnetic wave conformal to the 
upper interface of the 5-shape lens, although not exactly matches surface geometry especially 
for the realistic meta-atom made S"-shape lens, because of the inevitable scattering caused 
by the structure. Besides, it is worth noting that the electric field amplitude is not uniformly 
distributed in the space above the lens, because the electromagnetic wave will be focused 
by the concave part of the 5*-shape lens, but diffused by the convex part. Meanwhile, when 
the out-going electromagnetic wave just crosses the upper interface of the 5-shape lens, the 
simulation results of the average phase variation versus different frequencies [Fig. 5(c)] for 
5-shape lenses made of effective medium meta-atom (blue dots) and realistic meta-atom 
(black circle) indicate a near constant phase variation in the designed operating frequency 
range, implying the broadband ENZ response. And the mean value of the average phase 
variation is about 0.7rad (red line) in the operating frequency range. 



IV. CONCLUSIONS 

In conclusion, based on the Bergman spectral representation of the effective permittivity, 
we have proposed a design to achieve the realistic broadband ENZ meta-atom consisting 
of step-like metal-dielectric multilayered stacks in optical frequency range. The broadband 
ENZ response is achieved due to the fact that the electric field can penetrate through the 
meta-atom through different paths corresponding to different ENZ frequencies. Functional 
optical devices with exotic electromagnetic properties, including directional emission and 
phase front shaping, over a wide frequency range of more than 50 THz have been constructed 
with the meta-atom, and these unconventional devices can find many applications in optical 
communications, imaging processing, energy redirecting, and adaptive optics. Furthermore, 
our design can also be used in designing metamaterials with other broadband constant 
parameters, such as magnetic permeability, electrical conductivity, and thermal conductivity 
and diffusivity. 
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APPENDIX 

In order to achieve the best simulation results corresponding to the effective permittivity 
results depicted in Fig. 1 of the main text, the values of the gold inclusion filling ratio 
and the layer thickness for the realistic structures in Table 1 are calculated from the exact 
theoretical analysis and they are up to six significant digits mathematically. 

TABLE Al. Filling ratio of the gold inclusion and the thickness of each layer 



Layer 



fi 



di (nm) 



Theoretical Approximated Theoretical Approximated 

1 0.216392 0.22 43.6884 44 

2 0.188523 0.19 19.2691 19 

3 0.163955 0.16 37.0425 37 
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FIG. Al. (Color online) Robustness analysis of the designed parameters in the realistic meta-atom 
with solid curves indicating theoretical results in the main text based on six-digit values of the 
fraction and the thickness, while the dashed curves indicating the approximated results based on 
two-digit values. 
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In the practical fabrication, the feature size of the fabricated thin film has the accuracy of 
nanometer. Regarding the robustness of the current design, a simple analysis is performed 
based on the theoretical design, with the results indicated in the above Table Al and Fig. Al. 
As shown in Table Al, the fraction and the thickness values are approximated into only two 
significant digits. According to Fig. Al, the effective permittivity calculated with the ap- 
proximated two-digit values will only has a slight variation and small fluctuation, compared 
to the theoretical results obtained from the exact six-digit values in the manuscript. It 
is demonstrated that the performance of our design is quite robust to realistic fabrication 
tolerances. 

On the other hand, the impedance marching with free space is of importance in practi- 
cal applications. According to the obtained effective permittivity shown in Fig. 1(b), the 
impedance spectrum of the broadband ENZ metamaterial devices is plotted in Fig. A2. 
Within the broadband ENZ frequency region, the impedance also gives oscillatory response 
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FIG. A2. (Color online) The impedance of the designed broadband ENZ metamaterials 



with the definition of Z /Z = ^[i/e x , where Z is the impedance of free space. Although 
the impedance is not perfectly matched with free space, the value is relatively low within 
the ENZ frequency region, which is benefit from the existing of optical loss. 
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FIGURE CAPTIONS 

FIG. 1. (Color online) The schematic diagram of the broadband ENZ meta-atom and the 
designed broadband near-zero effective permittivity, (a) The effective-medium meta-atom 
with three homogenous layers composed by gold and silica, and the realistic meta-atom with 
a step-like metal-dielectric multilayer structure, (b) The theoretical values (solid curves) and 
the FDTD retrieved values (dashed curves) of the real part (red curves) and the imaginary 
part (blue curves) of the effective permittivity e x of the effective medium meta-atom, with 
the ENZ frequencies (yellow dots) listed at the right bottom. 

FIG. 2. (Color online) The distributions of the absolute value of the electric field in the 
effective medium meta-atom (upper panel) and the realistic meta-atom array (lower panel) 
at different ENZ frequencies, corresponding to the incident electromagnetic wave indicated 
by the coordinates. 

FIG. 3. (Color online) The schematic diagram and IFC analysis about the directional 
emission with respect to the broadband ENZ response, (a) The diagram of the directional 
emission through the ENZ prism represented by the electromagnetic waves inside and outside 
the prism, (b) The flat elliptical IFC (red curve for e x = 0.05 and e y = 3.0) of the meta- 
atom under lossless condition, (c) The elliptical-like IFC (red curve for e x = 0.05 + 0.5i and 
e y = 3.0 + 0.002i) of the meta-atom under lossy condition. The IFC of air is plotted as blue 
circle. 

FIG. 4. (Color online) Simulation results about the directional emission of the broadband 
ENZ prism with the I5°-oblique upper interface, represented by the distribution of the 
electric field amplitude E x and the power flow (hollowed arrows), for (a) the effective medium 
made prism and (b) the realistic meta-atom made prism, (c) The angle of refraction (in 
degrees) as a function of the frequency, based on the theoretical effective permittivity (blue 
curve) , and the simulation results of the effective medium made prism (blue dots) and the 
realistic meta-atom made prism (black circles). The average angle of refraction is 6.3° (red 
line) in the ENZ frequency range. 

FIG. 5. (Color online) Simulation results about the phase front shaping of the broadband 

12 



ENZ S-shape lens with the radius of curvature of 3 /zm, represented by the distribution of 
the electric field amplitude E x , for (a) the effective medium made S-shape lens and (b) the 
realistic meta-atom made S'-shape lens, (c) The average phase variation (in radians) as a 
function of the frequency, with respect to the effective medium made S'-shape lens (blue 
dots) and the realistic meta-atom made S-shape lens (black circles). The mean value of the 
average phase variation is 0.7rad (red line) in the ENZ frequency range. 



13 



/g\ Effective medium meta-atom 
■ d=100nm j 



!/(A'..l 



15 



2(H.) i(£„) 



*kH* 



Realistic meta-atom 
rf=100nn 



• :L ^^ H 


l 




SiO, 
3 


Au !*■ 


"Sr-i- 


P^ 1 - 


*i, 



*l K 



*KH- 



(b) 



0.4 

0.2 
0.0 

-0.2 
-0.4 
-0.6 



Imte) s 




/^ 




S^\ ii r/\i v 






i \/ iii \/ 


V 


: /rc(sj 


ENZ 


i. 537.42THz 
ii. 550.47THz 
iii. 563.52THz 
iv. 576.57THz 
v. 589.62THz . 



540 560 580 600 
Frequency (THz) 



FIG. 1. 



14 




y(K) 

537.421Hz 550.47TH; 

z (H„)<b > x (E ) 



FIG. 2. 



15 



(a) 



AIR 


e 


ft* 


ENZ 



z (H )cb > x (£.) 




fc 'M!: 



/ M 


i/0\ 




'\P*r^~- k " 


■' ~y- 


^v*- 


^^-^ \ ' 




■<~\ 


AIR^/_ 







FIG. 3. 



16 



yiK) 




537.42THz 



537.42THz 



x(EJ 



15 



10 



(c) 
Re(£j<0 


ENZ 

/~\ O O /\ 
/. \ /• \° t 


Re(£,)>0 

yv 


— -""■*" ° 

• o o 
■ . . .0 . . . . . 



> 





0L 

500 520 540 560 580 600 620 640 

Frequency (THz) 



FIG. 4. 



17 



R=3|jm 



s 



(b) 




( 



y(K) 



537.42THz 
2(H„) x(E ) 



537.42THz 



(c) 
Re(e,)<0 


ENZ 


Re(Ej>0 

® 


o ° o 


O 

• • 


O . . 

o 
) 





T3 

rj 

« 1-2 

o 

I 1.0 
g 0.8 
"I 0.6 

S3 0.4 

< 
0.2 

500 520 540 560 580 600 620 640 
Frequency (THz) 

FIG. 5. 



18 



